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Abstract 

The effect of inelastic final-state interactions (IFSI's) on the determination of the 
weak phase a from the isospin triangles of B — > tttt is qualitatively illustrated. Neglect- 
ing the electroweak penguins and IFSI's and assuming the dominance of the top-quark 
loop in strong penguin diagrams, we propose an experimentally accessible way to ap- 
proximately determine a and the penguin amplitude in — > tttt. This approach 
relies on a simplified isospin consideration and the factorization approximation, and 
its feasibility is irrelevant to the time-dependent measurements of Bj — > -kit. 
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1. Introduction 



Today the most promising way to test the Cabibbo-Kobayashi-Maskawa (CKM) mecha- 
nism of quark mixing and CP violation is to measure CP asymmetries in neutral 5-meson 
decays to CP eigenstates [p]]. For this purpose, the decay mode B d — * tt + tt~ is a good can- 
didate in addition to the gold-plated channel Bd — > J/ipKs- However, the penguin-induced 
effect on Bd — > tt + tt~ may be significant enough so that a reliable prediction for the CP 
asymmetry is theoretically difficult 0. Using the isospin triangle relations among the decay 
amplitudes of — > 7r ± 7r°, Bd ~ * 7r + 7r~ and 5^ ^ 7r°7r°, Gronau and London have shown 
that the tree-level and penguin contributions can be disentangled by measuring the relevant 
decay rates |§. Their work provides a relatively clean way to extract a, an angle of the CKM 
unitarity triangle in the complex plane || , although its feasibility depends closely upon the 
time-dependent measurements of Bd — > 7r7r. 

The naive isospin analysis of B — > tttt has to be modified, if these decay modes involve 
non-negligible electroweak penguin (EWP) diagrams H or inelastic final-state interactions 
(IFSI's) @. It has recently been shown that the overall effect of EWP's on the determina- 
tion of a from the isospin triangles is merely at the level of 0(X 2 ) (A = 0.22) or smaller 0, 
although the individual EWP amplitude may be 0(A) of the strong penguin amplitude in 
Bd — > 7T 7r°. The IFSI's are possible to induce channel mixing (e.g., between the direct reac- 
tions Bd — > tttt and the two-step processes Bd — > DD — > tttt), such that both the branching 
ratios and CP asymmetries of Bd — > tttt are difficult to be predicted. An illustration of the 
effect of IFSI's on the extraction of a from B — > tttt is still lacking in the literacture. 

In this work, we first illustrate the influence of tttt DD scattering on the isospin 
triangles of B — > tttt qualitatively. We find that the determination of a may suffer from 
large mixing between the 1 = states of tttt and DD. The approach proposed in ref. 
can work only if the IFSI's are insignificant enough to be negligible. Second, we develop 
an approximate but more practical approach to isolate the weak phase a and the penguin 
amplitude in B d — > tttt. Neglecting the effects of EWP's and IFSI's and assuming the 
dominance of the top-quark loop in strong penguin diagrams, we relate the tree-level and 
penguin amplitudes of B — ► tttt under a simplified isospin consideration and the factorization 
approximation. Provided that the decay rates of B^ — > tt^tt , B®^ tt + tt~ and B 2 ~~ > vr°7r° 
are measured, we show that it is possible to determine the magnitudes of the penguin and 
tree-level amplitudes and to extract the weak and strong phase shifts between them. This 
approximate method has two obvious advantages: (1) its feasibility depends only upon the 



time-independent measurements of the relevant decay rates, which can be carried out at 
either P-meson factories or high-luminosity hadron machines; and (2) it can confront the 
nearest data of B — > titt and give a ballpark number to be expected for a, before a delicate 
determination of a is available in experiments. 

2. Effect of IFSI's on the isospin triangles 

In the absence of EWP's and IFSI's, the isospin amplitudes of B d — > 7r + 7r~, P+ — > 7r + 7r°, 
B d — > 7T°7i° and their CP-conjugate counterparts are given as || 
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where some Clebsch-Gordan coefficients have been absorbed into the definitions of A^ \ and 
2- Clearly the above relations form two isospin triangles in the complex plane: 

A + _ + V2A 00 = V2A +0 , A + _ + V2A m = V2A_ . (2) 

Since the decay modes B^ —>■ 7i ± n occur solely through the tree-level quark diagrams, the 
ratio AY /AY = A-q/A+q is purely governed by the CKM factor (V ub V* d ) / (V* b V ud ) = e" 2i7 
(7 is an angle of the unitarity triangle ||). To probe CP violation induced by the interplay 
of decay and B d — B d mixing in B d — > 7r + 7r~, we need to measure the interference term 

Im£ + _ = Im l^- 2 ^ ^±=-\ = Im (e 2ia , (3) 

where f3 and a are also the angles of the CKM unitarity triangle (a + (3 + 7 = 180°), 
z = AY I AY and z = AY /AY- As shown in ref. 0, the complex parameter z (z) is 
determinable from the isospin triangles up to a two-fold ambiguity in the sign of its phase. 
In a similar way one can discuss the CP-violating term Im£ o i n B d — > 7r°7r°. Then the weak 
angle a can be extracted from Im£ + _ and Im£ o- This approach relies on the time-dependent 
measurements of B d ^ tttt. 

Now let us illustrate the effect of IFSI's on the isospin triangles and the CP-violating 
measurable Im£ + _ (Im£ o)> which are crucial for the extraction of a. Due to IFSI's the 1 = 



state of 7T7T can mix with that of DD, and this leads to the two-step decays Bd — > DD — > tttt. 
The final states 7r ± 7r° cannot mix with DD because of their different isospin configurations 
(1 = 2 for tt^tt , 1 = and 1 = 1 for DD). But in principle tt^tc could mix with p ± p°, 
and DD could mix with D*D* ||. For simplicity, here we only consider the influence of 
mi # DD scattering on Bd — ► tttt, leaving the 1 = 2 state of rnr unmixed with others. As a 
result, the "bare" isospin amplitudes of Bd — > tttt in eq. (1) are modified as 
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In the above equations, Aq D denote the 1 = amplitude component of Bd — > DD, 5*0,2 are 
complex matrix elements connecting the unitarized isospin amplitudes to the bare ones B. 
Of course, Sq D = and = Sq 7 " = 1 if there is no mixture between the 1 = states of 
mi and DD. From eq. (4), we obtain the following triangle relations: 

A' + _ + V2A' 00 = V2S™A +0 , A' + _ + V2A' 00 = V2S%*A- . (5) 

Comparing this result with eq. (2), one can observe an overall change of the bare isospin 
triangles due to IFSI's. In this case, the CP-violating measurable Im£ + _ becomes 

Im^_ = Im (e~ 2tf3 j^) = Im (e 2ia i— ^ , (6) 



where 
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We find that in principle the angle a can still be extracted from Im£ + _ (and Im£ 00 ). However, 
the determination of z' and z! from eq. (5) needs the knowledge of , since the measured 
branching ratio of B+ — > 7r + 7r° can only fix \A +0 \ other than |S , 2 T7r A +0 |. Before we obtain 
definite information on the mixing of tttt and DD, it is impossible to determine the values 
of \S^ D \ and \S^ 2 \ reliably []. Thus the isospin approach of ref. |§ cannot work well in the 
presence of significant IFSI's. 

If the 7i7x # DD scattering effect is insignificant in the problem under discussion, one 

expects that \Sq D \ may be small enough and jS'J 71 "! does not deviate too much from unity. In 

3 From ref. ||, one can obtain a rough feeling of the possible effect of IFSI's on the magnitudes of decay 
rates and CP asymmetries. 
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this case, the two triangles in eq. (5) are able to be approximately determined from direct 
measurements, then one can isolate the unknow parameters z' and z' to an acceptable degree 
of accuracy and extract a from Im£l_ and Ihi£q . In the following we shall assume this simple 
case and present an approximate approach to determine a and the penguin amplitude in the 
decay modes — > irir. 

3. Approximate isolation of the penguin amplitude 



To lowest order in weak interactions, all two-body mesonic B decays can be generally 
described by ten topologically different quark diagrams 0. In this language, we assume that 
the effect of IFSI's on every decay mode is insignificant enough to be negligible. Neglecting 
the EWP diagrams and those exchange- or annihilation-type channels, the dominant quark- 
diagram amplitudes for B® — > 7r + 7r~, B£ — > 7r + 7r° and B® — > tt°tt° are illustrated in fig. 1. 
Another reasonable assumption to be used is that the strong penguin diagram is dominated 



by the top-quark loop [JlOj, |TTJ. This implies that the weak phase of the penguin amplitude 
comes only from V^ h V t d in B® — > 7r + 7r~ or B® — > 7r°7r°. In terms of the angles of the CKM 
unitarity triangle ||, we parametrize the decay amplitudes of the three modes and their 
CP-conjugate counterparts as follows []: 
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The above expressions are based on a simple isospin consideration and the factorization 
approximation for the final states tttt, as we shall explain below. Here S denotes the strong 
phase shift between the penguin and tree-level amplitudes of B^ — > tttt. Thus without loss of 
any generality, T and P can be assumed to be real and positive. The parameter a represents 
the color- mismatched suppression in the tree- level amplitudes of B+ — > 7r + 7r°, B® — > 7r°7r° 
and their CP-conjugate processes. We shall see later that a can be either calculated using 



4 Here we use the same valence-quark notations for the pions as refs. jjJI, \n + ) = \ud), \tt°) = 

\dd — m)/V2, and 1 7r ) = — \du). This convention guarantees that the three pion mesons belong to an 
isotriplet. Also, we adopt \B° d ) = \db), \B° d ) = \bd), \B+) = \ub), and \B~) = -\bu). 



the factorization approximation or measured directly from the decay rates. Clearly the 
relations in eq. (8) can form the same isospin triangles as those in eq. (2). This implies 
that the quark-diagram language used for describing B — > -kit is consistent with the isospin 



analysis, as shown first in ref. [DJ. In the factorization approximation, we can assume that 
the difference between the interactions of / = and 1 = 2 final states is negligible. Thus the 
two quark-diagram amplitudes of — > 7T + 7r°, which occurs solely through the AJ = 3/2 
transition, can be parametrized as eq. (8a) with a real factor a. This accordingly implies 
that the tree-level amplitude of B® — > 7r + 7r~ and that of B® — > 7r°7r° also have negligible 
strong phase shift, to guarantee the isospin triangle relations in eq. (2). So do the penguin 
amplitudes of B® —>■ 7r + 7r~ and B® — ► 7r°7r°. Therefore, we expect that contributions from 
the penguin diagrams of Bd — > vr7r are crucial for the isospin triangles to have nonvanishing 
area. For a detailed discussion about translation of the isospin amplitudes of B — > tttt into 
the corresponding combinations of quark-diagram amplitudes, we refer the reader to ref. 
|T3. Based on the correspondence between the quark-diagram description and the isospin 



analysis, we subsequently use eq. (8) to relate the weak and strong phases in Bd — > 7T7t to 
the time-independent measurables. 

From eqs. (8a) and (8b), we have |A+.o| — |^4-o|- The decay rates of the above six pro- 
cesses should be measured in the near future at either 5-meson factories or high-luminosity 
hadron machines. Let us define the following four measurables: 



R +- = t~a — fo r 7 — r? ' R +- = ri — r?~m — I? > ( 9a ) 







12 


|A 


,o| 2 + 


A- \ 2 




l^oo 


2 


1^4- 


|2 _|_ 

r0| + 


A I 2 

^-o| 







2 




2 + 


A_ 


2 



R °° - U 12, 14 12 ' - \A 12 , M 12 • ( 9b ) 
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It is clear that T can be directly obtained from |4 +0 | 2 or |A_ | 2 . On the other hand, the 

magnitude of the penguin amplitude P is only determined once the data on (R+-) and 

-Roo (Roo) become available. By use of eqs. (8) and (9), the ratio P/T = \ is given as 

X = ^(1 + a) (ai? + _ + 2i? 00 ) - a = yj (1 + a) (aR + _ + 2R 00 ) - a . (10) 

It should be noted that the size of x is sensitive to the uncertainty associated with the color- 
suppression factor a. There are two ways to determine a. One way is to use the effective 
weak Hamiltonian 7i(AB = 1) |14| and the factorization approximation [15]. It is easy to 
obtain 

where iV c is the number of colors, C% ~ —0.291 and Ci ~ 1.133 are the Wilson coefficients 
at the scale \i = 0(m b ) @. We obtain a w 0.23 for A^ c w 2.2 and a w -0.26 for N c = oo. 



The current data on exclusive hadronic B decays favor the former, i.e. a positive a with 
the magnitude of 0(A) ||16|| . The other way to determine a is through measurements of the 
decay modes Bd — > ti + ti^ and Bd — > 7i tt°. From eqs. (8) and (9) one can find 

a = -2- -= — . (12) 

rt-j — it^ 

Indeed the rate differences _R o — -Roo an d R+- — R+- indicate the signals of direct CP 
violation. The assumptions made in eq. (8) imply that the direct CP asymmetries in 
Bd — > 7r + 7r - and Bd — > 7T 7r° are governed by the same weak phase difference {(3 + 7) and 
the same strong phase shift (5), but their magnitudes are different. A comparison between 
the results of a obtained in eqs. (11) and (12) can serve as a rough test of the factorization 
approximation. In practice, however, measuring the rate difference -Roo — Roo might be very 
difficult due to the expected smallness of the branching ratio of Bd — > tt tt . Considering 
SU(3) symmetry and its small breaking effect, Deshpande and He |T7| have shown that the 



rate difference of B® — > 7r + 7r~ (7r°7r°) vs B® — > 7r + 7r~ (7r 7r ) is related to that of B^ — > n~K + 
(n K ) vs B® tt + K~ (tt°K°). The latter can be more easily measured in experiments, 
thus can be used to determine a through eq. (12). It should be noted that our main results 
still hold even in the special case a = 0. If S happens to be 0° or 180°, however, eq. (12) will 
not work. In this case, we have to use the data on other B decays to determine a fT6|j . 

4. Approximate determination of a 

Now we look at how to extract the weak phase a without resorting to measurements of 
the decay-time distribution of Bd — ► 7r7r. Since the weak phase shift between the penguin and 
tree-level amplitudes of Bd — > is /3-py, by use of the unitarity condition a+/?+7 = 180° 

we obtain two concise equations connecting the weak and strong phases to the measurables: 

cos(a + 5) = — \l + x 2 ~ (1 + a) 2 R + J\ , (13a) 
2x L J 

cos(a -5) = — [l + x 2 ~ (1 + a) 2 R + ] . (13b) 
2% L J 

Similarly, a and 5 can be related to i?oo and -Roo- Because x is determined from eq. (10), 

both a and 5 can be extracted from (13) with four-fold ambiguity. Due to the assumptions 

made in eq. (8), this ambiguity cannot be reduced by the relations among a, 5, Rqq and 

-Roo- To the accuracy that and _R + _ cannot be distinguished in experiments, we expect 

that the strong phase shift 5 might be close to 0° or 180° (i.e., vanishingly small direct CP 

violation in Bd — > titt). Note that eq. (13) would become invalid if the magnitude of the 
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penguin amplitude P were too small to be measured. In this case, the two isospin triangles 
in eq. (2) would collapse to lines. Within the standard model, however, several calculations 
have given that x is a t the level of 0(A) or larger 0. Thus our results (10) and (13) are 
feasible to probe the weak phase a as well as the significant penguin amplitude in B d — > nn. 



It is worth remarking that our approach to determine the angle a does not require the 
measurement of the CP- violating term Im£ + _: 



sin(2a) — 2% sin a cos S 



(14) 



(to lowest-order approximation of \) > which is induced by the interplay of decay and B d — B® 
mixing. Thus the feasibility of our approach is independent of the time-dependent measure- 
ments of neutral B decays to be carried out at the asymmetric B factories of KEK and SLAC 
r8|. Nonetheless, it should be noted that the determination of and \A\ | from any 



experiment is indeed influenced by B® — B® mixing. To overcome this problem, one available 
way is through the time- integrated measurements of Bd — > 7r + 7r~ on the T(4S*) resonance, 
where the produced two Bd mesons are in a coherent state (with odd charge-conjugation 
parity) until one of them decays. Using the semileptonic decay of one Bd meson to tag the 
flavor of the other meson decaying to 7t + tt~, the probability for observing such a joint decay 
has been give as |T9|: 



Pr{rX + ;n + n-) = 
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|(Z X + \H\B d v )\, and Xd = Am/T w 0.71 is a measure of 



B® — B^ mixin. 



l + X~ and B Q d 
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J]. At present, the semileptonic decays B® 
have been well reconstructed || 2~0| , thus \ A{\ is determinable independent of the above joint 
decays. Once Pr(Z ± X =F ; 7r + 7r _ ) are measured, we shall be able to determine \A + _ \ and 
through eq. (15). In practice, the measurements of Pr(Z ± X =F ; 7r + 7r _ ) can be realized either 



by the symmetric e + e collider running at Cornell |2(| or by the forthcoming asymmetric B 
factories. Similarly, \A QQ \ and \A m \ can be determined. 



5. Discussions and conclusion 



In the first part of this work, we have qualitatively illustrated the effect of IFSI's on the 
isospin triangles of B — > tttt. The determination of a from these decay modes may suffer 



s 



from large mixing between the 1 = states of tttt and DD. A quantitative estimation of 
IFSI's needs much progress in theory and more data in experiments. 

In the assumption that the IFSI's in B — > titt are insignificant enough to be negligible, 
we have presented a self-consistent and time-independent approach to approximately isolate 
the weak phase a and the penguin amplitude of B<j — > 7T7T. Here it is worthwhile to comment 
briefly on the uncertainties with this approach. In the first approximation, we neglected 
the electroweak penguins in B — ► titt. As pointed out in ref. J7|, the effects of electroweak 
penguins on the decays via b — > uud or b — > d are at most 0(A 2 ) of the dominant tree-level 
amplitude T or 0(A) of the strong penguin amplitude P. Hence it is unlikely that the validity 
of the isospin analysis in ref. |J and our present method are significantly affected by the 
neglect of the electroweak penguins. On the other hand, the approximation of neglecting 
those exchange- and annihilation-type topologies in B — > is expected to be safe for some 



dynamical reasons fL3 |. The approach to test the reliability of this kind of approximations 



has been suggested in ref. The factorization approximation used to parametrize the 

decay amplitudes of B — > nn as eq. (8) is equivalent to neglecting the difference between 
the 1 = and 1 = 2 interactions in the final states. The validity of this assumption can be 
examined by the forthcoming experimental data. The assumption that the strong penguin 
diagram is dominated by the top-quark loop can be interpreted as follows. On the basis of 
the standard model with three families of quarks, the penguin contributions with u and c 
quarks in the loop are equal in magnitude up to 0(m 2 c /rn%) apart from their CKM factors. 
Thus the top-quark loop dominates the penguin amplitude to the same degree of accuracy, 
i.e. 0(< 10%), guaranteed by the CKM unitarity. Indeed this is a good approximation to 
obtain the lowest order magnitudes of a, 5 and P. When the time-dependent measurements 
of Bd — > hit become available at an asymmetric B factory to allow for the more precise 
isospin analysis proposed by Gronau and London ||, one can check the difference between 
our results and theirs. The corresponding deviation should definitely constrain the reliability 
of the top-quark dominance and the factorization approximation assumed in our work. 

By now the experimental sensitivity to the penguin contribution in B^ — > hit has been 
analyzed based on the isospin triangle relations || and under the circumstance of asymmet- 
ric B- meson factories [|ST]. Such an analysis is also applicable to the simplified approach 



proposed in this work. It is concluded that determining the weak and strong phases associ- 
ated with the penguin amplitude to a satisfactory accuracy is experimentally feasible in the 
near future. Finally we stress again that our approach is time-independent, and thus it can 
be confronted with all the forthcoming measurements of B — > tttt, no matter whether they 



o 



are time- integrated or time-dependent. 



Acknowledgements 

We are very grateful to X.G. He, T. Nakada and S. Stone for useful discussions and 
comments. Z.Z.X. would like to thank H. Fritzsch for his warm hospitality and constant 
encouragement. C.H. Would like to thank N. Khuri and the Physics Department of Rock- 
efeller University, where part of this work was done, for their warm hospitality. The work 
of C.H. was partially funded by the N.S.E.R.C. of Canada, and the research of Z.Z.X. was 
supported by the A.v.H. Foundation of Germany. 



References 

[1] A.B. Carter and A.I. Sanda, Phys. Rev. Lett. 45 (1980) 952; 

LI. Bigi and A.I. Sanda, Nucl. Phys. B193 (1981) 85; B281 (1987) 41; 

for a review, see: I.I. Bigi, V.A. Khoze, N.G. Uraltsev, and A.I. Sanda, in CP Violation, 

edited by C. Jarlskog (World Scientific, Singapore, 1989), p. 175. 

[2] M. Gronau, Phys. Rev. Lett. 63 (1989) 1451; Phys. Lett. B300 (1993) 163; 
D. London and R. Peccei, Phys. Lett. B223 (1989) 257; 
B. Grinstein, Phys. Lett. B229 (1989) 280; 
D. Du and Z.Z. Xing, Phys. Lett. B280 (1992) 292. 

[3] M. Gronau and D. London, Phys. Rev. Lett. 65 (1990) 3381; 

H.J. Lipkin, Y. Nir, H.R. Quinn, and A.E. Snyder, Phys. Rev. D44 (1991) 1454. 

[4] C. Hamzaoui, J.L. Rosner, and A.I. Sanda, in the Proceedings of the Workshop on 
High Sensitivity Beauty Physics at Fermilab, edited by A.J. Slaughter et al. (Fermilab, 
Batavia, IL, 1988), p. 165. 

[5] N.G. Deshpande and X.G. He, Phys. Rev. Lett. 74 (1995) 26; preprint OITS-572 (1995). 

[6] M. Wanninger and L.M. Sehgal, Z. Phys. C50 (1991) 47; 
A.N. Kamal, Int. J. Mod. Phys. A7 (1992) 3515. 

[7] M. Gronau, O.F. Hernandez, D. London, and J.L. Rosner, preprint |hep-ph/9504327 . 

[8] Particle Data Group, L. Montanet et al., Phys. Rev. D50 (1994) 1173. 

1 n 



[9] Z.Z. Xing, preprints [hep-ph/9502339| and [hep-ph/9409318| (both to appear in IL Nuovo 
Cimento A). 

[10] The same assumption has been used by several authors. See, e.g., 

C. Hamzaoui, A.I. Sanda, and A. Soni, Phys. Rev. Lett. 63 (1989) 128; 

Y. Nir and H.R. Quinn, Phys. Rev. Lett. 67 (1991) 541; 

J. P. Silva and L. Wolfenstein, Phys. Rev. D49 (1994) R1151. 

[11] M. Gronau, J.L. Rosner, and D. London, Phys. Rev. Lett. 73 (1994) 21. 

[12] D. Zeppenfeld, Z. Phys. C8 (1981) 77. 

[13] M. Gronau, O.F. Hernandez, D. London, and J.L. Rosner, Phys. Rev. D50 (1994) 4529. 
[14] A.J. Buras et al, Nucl. Phys. B370 (1992) 69. 

[15] A.J. Buras, J.M. Gerard, and R. Riickl, Nucl. Phys. B268 (1986) 16; 
B. Blok and M. Shifman, Nucl. Phys. B389 (1993) 534. 

[16] M.S. Alam et al., preprint CLNS 94-1270 & CLEO 94-5 (1994); 

T.E. Browder, K. Honscheid, and S. Playfer, in B Decays (2nd edition), edited by S. 
Stone (World Scientific, 1994). 

[17] N.G. Deshpande and X.G. He, preprint OITS-566 (1995). 

[18] T. Nakada, preprints PSI-PR-91-24 and PSI-PR-92-11; 
Y. Sakai, Nucl. Inst. & Meth. A333 (1993) 42; 
BaBar Collaboration, preprint SLAC-443 (1994); 
Z.Z. Xing, preprint |hep-ex/9502001| (to appear in Phys. Rev. D). 



[19] D. Du and Z.Z. Xing; Phys. Rev. D47 (1993) 2825; 

H. Fritzsch, D.D. Wu, and Z.Z. Xing, Phys. Lett. 328 (1994) 477; 
Z.Z. Xing, Phys. Rev. D50 (1994) R2957. 

[20] S. Stone, Nucl. Inst. & Meth. A333 (1993) 15; 

S. Playfer and S. Stone, preprint HEPSY 95-01 (to appear in Int. J. Mod. Phys. Lett.). 

[21] R. Aleksan, A. Gaidot, and G. Vasseur, in the Proceedings of B Factories — the State of 
the Art in Accelerators, Detectors and Physics, edited by D. Hitlin, SLAC-400 (1992). 



1 1 




u 
u 



TT 



+ 



TT 



u - 




u 
d 



TT 



7T 



+ 



(bl) 



(b2) 




(cl) 



u 
d 



TT 




(c2) 



TT 



TT 



Figure 1: The dominant tree-level and strong penguin diagrams for — > tt + tt , 

7T + 7T , and B° d -> 7T 7T . 
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